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Abstract 

Suppose A(x), B(x) are 2x2 matrices on an interval [0, oo) and C a constant diagonal matrix 
with distinct positive entries. Let U(x,t) be the matrix solution of the system of hyperbolic 
PDEs CUu - Vxx - AU X - BU = on [0,oo) x R with the initial condition U(-,t) = for 
t < and the boundary condition U(0,t) = 6(t)l2- We prove a stability result for the inverse 
problem of recovering A, B from U x (0, •). The solutions of the forward problem propagate with 
two different speeds so techniques for inverse problems for a single hyperbolic PDE are not 
applicable in any obvious way. 



1 Introduction 



Below, =<! will mean an inequality up to a constant multiple, all functions will be real valued, upper 
case letters such as M will represent 2x2 matrices with entries Mij, lower case bold letters such as 
v will represent 2x1 vectors with components vi , x>2 ■ All convolutions will be in the t variable if the 
convolution involves a function of x and t. We define the operator C by Cv := Cvu— v xx — Av x — Bv 



where C 



A 2 




with < A < fi, A(x),B(x) are real valued 2x2 matrices and v(x,t) is a 
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2x1 vector. Let U(x,t) be the real valued matrix solution of the IBVP 

CU = for (x, t) G [0, oo) x R (1.1) 
U = for t < (1.2) 
J7(a; = 0, t) = 5{t)I 2 for t G R. (1.3) 

where -Z2 is the 2x2 identity matrix. We study the recovery of A(-),B(-) or a subset of these 
coefficients if we are given U x (x = 0, •). 

Such an inverse problem arises in the examination of the structural integrity of a composite 
beam; please see the introduction of [MNS05] for a discussion of this application and also for 
other references related to this application. The problem of determining the (spatially varying) 
parameters for the Timoshenko model of a beam (see [A73]), from measurements of the deflection 
from the neutral axis and the twist in the cross-section, also may be modeled as an inverse problem 
for a two speed second order hyperbolic system on an interval with two dependent variables; the 
entries of A, B are made up of the parameters in the Timoshenko beam model. 

For future use we define u, u to be the columns of U, that is U = [u, u]. Then u and u also 
satisfy (|l.ip . (jl.2p but satisfy the boundary condition 

u(x = 0, t) = 5(t)ei, u(x = 0,t) = 5(t)e 2 (1.4) 

where ei and e 2 are the columns of I 2 . There are two speeds of propagation associated with 
C, namely 1/A and 1/fi and U\,U\ are the fast moving components and u 2 ,u 2 the slow moving 
components of u and u respectively. It is this feature of the problem which makes it difficult to 
apply any obvious modification of the inversion schemes popular for inverse problems for a single 
hyperbolic PDE in one space dimension. 

In Theorem[3]we show that fjl . f)l .3[) has a unique solution in C 2 ([0, oo),P'(R)) and we give a 
progressing wave expansion of U. We postpone the statement of the theorem about the existence 
and the structure of U to the end of this section since the statement is quite long and follows 
from the standard progressing wave expansion technique; we want to draw attention to the more 
interesting results in Theorems [1] and [2] stated below. 

Let D = diag(A) be the diagonal matrix formed by taking just the diagonal entries of A. Define 
the diagonal matrix M(x) := e~2 Jo D (y) d v an d define v := M _1 v. Then we may show that 

Cv tt - v xx - Av x - Bv = M(Cv a ~ v xx - Av x - Bv) 

where A = M~ l (A-D)M and B = M~ 1 (D 2 /4 — D'/2-AD + B)M. Note that the diagonal entries 
of A are zero. Further, M(0) = I so 17(0, t) = U(0,t) = 6(t)I and {M~ l U) x (0, t) = U x (0,t). So 
for every pair (A, B) one can construct a pair (A, B) with the same data (M~ 1 U) x (0, t) as (A, B) 
except that the diagonal entries of A are zero. Hence, below we will study only the situation where 
the diagonal entries of A are known. 

Define the operator C T by 

£ T v := Cv tt - v xx + (^v), - B T v = Cv tt - v xx + A T \ X - (B - A') T v. 
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If v(x, t) and w(x, t) are 2x1 vectors then one may show that 

v T £w - (£ T v) T w =(v T Cw t - vJCw)i + (vjw - v T w x - v T 4w), (1.5) 

implying C T is the formal adjoint of C. Hence C will be formally self-adjoint iff A T = —A and 
B T = B — A', that is iff the diagonal entries of A are zero and B — B T = A'. 

An analysis of the linearized inverse problem with the linearization done around A = 0, B = 
gives an indication of the results one may expect for the inverse problem under consideration. When 
A = 0, B = the solution of (UJ- CCl is 



U(x,t) 



S(t - Ax) 
S(t - nx) 



;i.6) 



Hence the linearized forward problem about the trivial background is the solution of the IBVP 

C(5U) tt - (5U) XX = (5A)U X + (SB)U, (x, t) e [0, oo) x R (1.7) 
(SU)(0,t) =0, SU = for t < 0. (1.8) 

Here we assume that diag(5A) = 0. 

Fix a t > 0. We use (ll.5p with C corresponding to A = 0, B = 0, v(x, i) = U(x, r — t) and 
w(x,t) = (5C/)(x,t). Integrating this relation over the region [0, oo) x R, integrating by parts and 
using (|1.8j) . (|1.6j) . we obtain 

/oo 
[/(0,r-t) T (5C/) a; (0,t)^. (1.9) 

Now, using (|1.6p in (jl.9p and integrating one may show that 



(6U x ) n (Q,T) = ±(8B)iiMt)), (8U x ) 12 (0,t) = -L-(SB) 12 (x m (r)) - _-^-_(^)' la (x m (r)) 



(<5C4) 21 (0,t) = x l_(«5B) 2 i( a : m (r)) - ^ ^(6A)' 21 (x m (r)), (5U x ) 22 (0,r) = ±-(6B) 22 (x s (t)) 
where 

are the lengths probed from the origin, in time r, by a round trip using two fast waves, a fast and 
a slow wave, and two slow waves respectively. 

So, for this linearized inverse problem, where (5U) X (0, t) is given and 5 A, SB are to be deter- 
mined, one recovers the combinations SBn, SB^ — j^dA'^, 5B 2 i — j^5A' 2 i and 5B 22 . Hence if one 
is given two linearly independent relations amongst 5B\ 2 , 5B 2 i, 5A' 12 , SA 2 i, which are independent 
of 5B\ 2 — j^dA'^ = and 5B 2 \ — j^5A' 21 = 0, then one can recover SA, 5B from (5U) X (0, ■). For 
example, if we are given the value of 5 A then one can recover 5B. When the system is self-adjoint 
we have SB — (SB) T = SA' , that is SB\ 2 — SB 2 \ = SA' l2 = — SA' 2l . However, these relations are not 
independent of the two relations mentioned above, so we need an additional relation or the value 
of one of SBi 2 ,SB 2 i,SA' 12 , SA' 2l would have to be part of the data given. 
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Fi gure 1: Definition of x s , x m , xj 



This analysis suggests that, for the original inverse problem, given U x (0,t) on an interval [0, r] 
and the diagonal entries of A(-), one may expect to recover only four out of the remaining six 
coefficients in A, B, provided the other two coefficients are given. Further, the values of these 
coefficients will be recovered over intervals of different lengths which suggest that there may be 
complications using the downward continuation method popular for inverse problems for a single 
hyperbolic PDE in one space dimension. However, if all the coefficients except 622 are known then 
there should be no difficulty recovering 622 with the use of a downward continuation method. 

Our main result is a stability result for the original inverse problem and the proof reflects the 
discussion above. An examination of the proof will show that one may prove stability in more 
situations than covered in the statement of the theorem. 

Theorem 1 (Stability). Fix positive constants X andK. Suppose A, A G C 2 [0,X], B,B G C l [0,X] 
with \\A\\ C 2 + \\B\\ c i < K, \\A\\ C 2 + \\B\\ c i < K, and 

A(0) = 1(0), diag(A(-)) = = diag(A(-)). 

Let U and U be the solutions of U. 1\) - [T73\) corresponding to A,B and A,B respectively, on the 
region {(x,t) : < x, t + Xx < XX}. If either A(-) = A(-) or the off-diagonal entries of B(-) and 
B(-) are the same, then 

\(B-B)(x)\ + \(A' -A')(x)\ 4 max 1^(0, •) -U x (0,-)\, Vx G [0, XX/fj] (1.10) 

[0,2/iis] 

with the constant determined only by X , K , X and (i. 

The theorem suggests that given U x (0, t) over the interval [0, T] one should be able to reconstruct 
(some of) the coefficients over an interval [0, T/(2/i)] - the interval determined by the slower speed 
of transmission. Using the ideas discussed earlier, one may derive a result similar to Theorem Q] if 
the hypothesis diag(vl) = = diag(^4) is replaced by the weaker hypothesis diag(A) = diag(^4). 

For a p(t) G C 2 (M) with support in [0,oo), let v(x,t) be the solution of the IBVP 

£v = in [0, 00) x R 
v(0, t) = p(i) for t G R, v(-, t) = for t < 0. 
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The fastest speed of propagation being 1/A, it is clear that v(x,t) will be supported in the region 
< Xx < t. However, for certain choices of p(-), due to cancelations, the support of v(x,t) may 
lie in the slow region < fix < t. In [BBI97] Belishev et al made an important discovery where 
they showed that, if C is formally self-adjoint, then there is a unique function l(t) (independent of 
p(i)) so that if p\ = I * P2 then v(x,t) is supported in the slow region < fix < t. In fact, since 
v = pi * u + p2 * u (with the convolution in t alone), l(t) is the unique function so that u + / * u 
is supported in the slow region < fix < t. Using some of the ideas in [BBI97], we have extended 
their result to the general A, B case and simplified the proof. 

Theorem 2 (Existence of slow waves). If A 6 C 2 [0, oo) and B G C^fO, oo) then there exists a unique 
/(•) in C 2 [0,oo) so that u(x,t) + l(t) * u(x,t) is supported in the region < fix < t. Further, for 
any r > 0, ||/||c[o,r(^-A)/(/t+A)] bounded by a constant determined only by A, p and ||^4.||c[0,a:f(r)]> 



In [BBI97] and [BI02], Belishev et al studied the inverse problem considered in this article (for 
smooth coefficients though their arguments are valid for less regular coefficients) except with the 
additional requirement that C be formally self-adjoint. In this case there are only four coefficients 
to be determined but then U x (0, •) is also symmetric in this cas^ll so the data U x (0, •) consists of 
only three functions. With this in mind, Belishev et al in [BBI97] and [BI02], for the self-adjoint 
case, studied the recovery of B from U x (0, •) and /(•). They showed that B (and hence ^4) could be 
reconstructed from U x (0, •) and /(•). Further (for the self-adjoint case) they characterized the range 
of the map (A(-), -£?(•)) *— > U x (0, •); they showed that a function r{t) is in the range of this map iff a 
certain integral operator, defined in terms of r(t), is positive definite. Their proof showed that any 
pair of functions (r(t),l(t)) defined over appropriate intervals, with r(t) satisfying the "positivity 
property" is generated, in the above sense, by some A(-),B(-) associated with a self-adjoint L. 

Since l(-) is not an experimentally measurable quantity, in [BI03], again for the self-adjoint case, 
and assuming A was known, Belishev et al studied the recovery of B (three unknown quantities) 
from U x (0, •). They showed that they could reconstruct l(t), at least over a small interval, and hence 
from [BBI97] they could recover B over a small interval. Using this result Morassi et al in [MNS05] 
showed that if A = and B is symmetric (part of self-adjoint case) then the map B i— > U x (0, ■) 
is injective (uniqueness in the inverse problem). Our Theorem 1 covers the uniqueness (but not 
the reconstruction) results in the above references and we provide a fairly simple proof of stability 
for a more general situation. Belishev et al use the Boundary Control Method which has proved 
effective for reconstructions for several inverse problems for hyperbolic PDEs and Morassi et al 
combine this with a downward continuation argument in the frequency domain. We do not have 
a reconstruction method even if /(•) is part of the data. Finally, [Ni91] is a good starting point to 
read about the results of Nizhnik and his school on inverse problems for two velocity systems. 

Our proof of Theorem Q] uses a trick similar to the one used to analyze the linearized inverse 

For any r > 0, using (|1.5|l with w(s, t) = u(a;, t) and v(a;, t) = u(a;, r — t) and integrating over [0, oo) x R (using 
w(0, t) = 5(t)ei, v(0, i) — S(t — t)e2) we obtain 



B \\c[0,x f (r)]- 



= 




eIu*(0,t)(5(T -t) - uj(0,r - t)ei5(t) +e%A ei S(t - r)S(t) dt = d x Ui 2 (0,r) - d x U 2 i(0,T). 
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problem above. This trick was first used (as far as we know) in [SnSy88] for a single hyperbolic 
PDE and then applied to a system of hyperbolic PDEs in [Sa86], [SaSy87]. 



The existence and uniqueness of a weak solution of (|1.X[> - (|1.3|) may be proved by appealing to 

standard results but proving higher order piece-wise regularity requires dealing with some quirks 
in two speed problems. The following proposition characterizes the principal singularities in u and 
u and the existence theory associated with this expansion. 

Theorem 3 (Well posedness of the forward problem). If A £ C 2 [0, oo), diag(A(-)) = and 

B G C^C^oo), then there exist unique solutions u(x,t), u(x,t) in C 2 ([0, oo), V'(R)) of ( fi.ffj) . 
jXg ). Further, for all (x,t) £ [0, oo) x E, 

u(x, t) = S(t - Ax)ei + f (x, t)(H(t - Ax) - H(t - fix)) + g(x, t)H{t - fix), (1.11) 

u(x, t) = S(t - fix)e 2 + f (x, t)(H(t - Ax) - H(t - fix)) + g(x, t)H(t - fix), (1.12) 

where f, g are 2x1 column vectors which are C 2 solutions of the characteristic IB VP (see Figure 

w 

£f = in < Ax < t < fix, Cg = in < fix < t, (1.13) 
with the boundary, characteristic and transmission conditions 

g(0,t) = 0, t>0, (1.14) 

1 f x A f x 

fi(x,\x) = —J bn(z)dz+ 2 _ ^ J ai 2 {z)a 2 i{z) dz, (1.15) 

f 2 (x, Ax) = ^ 2 ^ 2 a 2 i(x), ((A 2 + fL 2 )f 2t + 2Xf 2x + Aa 2 i/i)(x, Ax) = b 21 (x), (1.16) 

(5i - fl)(x, fix) = 0, (gi - fi)t(x, fix) = p^2pa2i(0)ai 2 (x), (1.17) 

(92 ~ f2)(x,fix) = -^A_ a21 (o). (1.18) 

Further f, g are C 2 solutions of the characteristic IB VP (see Figure^) 

CA = m < Ax < t < fix, Cg = in < fix<t, (1.19) 
with the boundary, characteristic and transmission conditions 

g(0,t)=0, t>0 (1.20) 

/ 1 (x J Ax) = ^^ai 2 (0), (1.21) 

f 2 (x, Ax) = 0, f 2t (x,Xx) = jj2~^j^2a 12 (0)a 21 (x) (1.22) 
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Figure 2: Domains of f , f , g, g 



(91 ~ fi){x,fix) = 2 _ 2 ai 2 (x), 



(1.23) 



(92 ~ / 2 )(x,/xx) = a 12 (z)a 2 i(z) dz + — / b 2 2(z) dz 




(1.24) 



((A 2 + n 2 ){g\ - h) t + 2fi(gi - f x ) x + fia 12 (g 2 - f 2 ))(x,(j,x) = b 12 (x). 



(1.25) 



Using the ideas discussed earlier, one may derive a result similar to Theorem [3] if the hypothesis 
diag(A) = is dropped. 

It would be reasonable to ask if results similar to Theorems 1,2,3 hold if the boundary condition 
(jl.3p is replaced by U x (0, t) = 5(t)I 2 and for Theorem [1] the data is U(0,t) instead of U X (Q, t). We 
see no reason why the same methods will not work after adjusting the order of the singularity 
in U(x,t), that is the most singular term in the expansion of U(x,t) would be H{t — \x) and 
H (t — /ix) instead of 5(t — Xx) and 5(t — fix). 

The rest of the paper consists of the following. In section [2] we prove Theorem [TJ In section 
[3] we prove Theorem (2) Our proof uses some of the ideas in [BBI97] for the self-adjoint case, but 
we do not use the Boundary Control Method machinery and we think perhaps our proof is more 
transparent. In section H] we prove Theorem [3] and Proposition U] which is needed to complete the 
proof of Theorem [3j The proof of Theorem [3] consists of two parts : a progressing wave expansion 
and a well-posedness theory for a characteristic transmission boundary value problem for a system 
of equations. The progressing wave expansion part is standard but since the expressions are not 
in the literature we give the expressions and the derivation. The well-posedness theory for the 
characteristic transmission boundary value problem for a system with two velocities is not given in 
the literature though its proof uses standard techniques except for the appearance of an unusual 
transmission BVP problem for a single hyperbolic pde. 

Finally we wish to thank Mikhail Belishev for discussions about the problem considered in this 
article. 
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2 Proof of Theorem [T] 



Extend A, A as C 2 functions and B,B as C 1 functions, on [0, oo), with compact support, so that 
the C 2 norms of A, A and the C 1 norms of B,B, on [0, oo), are bounded by a constant multiple 
of the corresponding norms on [0, X], with the constant independent of A,A,B,B. Let U = [u, u] 
and U = [u, u] be the solutions of (|l.ip - (jl.3p corresponding to A, B and A, B respectively, over the 
region [0, oo) x R guaranteed by Theorem [3j Further, let l(-) and /(•) be the functions guaranteed 
by Theorem [2] for the operators corresponding to A, B and A,B. Note that the value of U x (0, t) 
and U x (0, t) for t £ [0, 2XX] is not affected by the extensions of A, A, B, B because the fastest speed 
of propagation is 1/A. 

Define 5A := A — A, 5B := B — B, 5U := U — U, (8a)ij := a^- — o^-, (5b) ij := b{j — bij, and 
\M\ := maxjj \rriij\. Note that the diagonal entries of 5 A are zero because of the hypothesis. We 
will prove the stability by showing a Volterra type estimate 

\(5B--^—(5A)'VC)(x)\ 
A + /x 

4 max \(5U) X (0, -)|+ / \(5U) x (0,t)\dt + \(5A)'(y)\ + \(5B)(y)\dy, x G [0, XX/ /i] 

with the constant determined only by A, /i, X and K. Then Theorem [1] follows from GronwalFs 
inequality and the hypothesis that either 5 A = or the off-diagonal entries of 5B are zero and the 
diagonal entries of 5 A are zero. 

The progressing wave expansions of u, u are given by (jl.lip . (|1.12p and from Theorem [3] 

u(x,t) = 5(t - Ax)ei +f(H(t - Xx) - H(t - /ix)) + gH(t - /ix) (2.1) 
u(x, t) = 5(t - fix)e 2 + f (H(t - Xx) - H(t - fj,x)) + gH (t - fix) (2.2) 

with f and g having properties similar to f , g. From Theorem [21 Theorem [3] and Proposition [J] we 
have that the C° norm of /, / on any finite interval and the C 2 norms of f , f , g, g on appropriate finite 
regions will be bounded by functions of A, //, X and K and parameters determining the interval or 
the region. Since the regions of interest below will be determined by A, \x and X, one is assured 
that all these norms are bounded by functions of A, /i, X and K. 

We will use the following four pairs of vector functions a(x, t), f3(x, t) defined on [0, oo) xK- 
I. a(x, t) = u(x, t) — u(x, t) and /3(x, t) = u(x, t — t); 

II. a(x, t) = (u — u)(x, t) + l(t) * (u — u)(x, t) and f3(x, t) = u(x, r — t) + (/ * u)(x, r — t); 

III. a(x, t) = u(x, t) — u(x, t) and /3(x, t) = u(x, r — t) + (I * u)(x, r — t); 

IV. a(x, t) = (u — u)(x, t) + l(t) * (u — u)(x, t) and f3(x, t) = u(x, r — t). 



For each of these pairs we note that 



S 



a(0, t) = on R and a(x, t) = for t < 

jCyS = on [0, oo)xl and 0(-,t) = for t » 0. 



Hence using (|1.5|) we have 

/'OO /"CXO /"CXO /'OO /'OO 

/ / f3 T £adtdx= / / (3 T £a - {£(3) T adtdx = / 0(0, t^a^O, i) eft. (2.3) 

j0 J —oo JO J — oo J —oo 



In each of the four cases 



a^O, t) = (u - 11^(0, t) or (u - 5)^(0, t) + / * (u - u) x (0, t), 
/3(0, i) = 5(t - t)ei or S(t - t)e 2 + Z(r - i)ei. 



Hence 



| RHS of | ^ \(U -U) x {0,t)\ + [ \(U-U) x (0,t)\dt (2.4) 

JO 

with the constant determined only by X and if. 



Estimating the LHS of (|2.3p . in each of the cases, may involve one of the following estimates for 
2x1 vectors v(x,t),w(x,t) which are C 1 and a continuous 2x2 matrix M[x). The derivation of 
these estimates is fairly straightforward with an integration by parts required for the first estimate. 

poo poo p — ^ — 

/ / v T Mw5'(T-t-a 1 x)H(t-a 2 x)dtdx 4 \M( )\+ T 1 ** \M(x)\dx; (2.5) 

Jo J-oo (Ji+a 2 Jo 

/•oo poo 

/ / v t Mw5(t -t- <rix)5(t- a 2 x)dtdx ^ |M( )|; (2.6) 

Jo J-oo <T\+a 2 

poo poo p — ^ — 

/ / v T Mw5(T-t-a lX )H(t-a 2 x)dtdx 4 j "* \M(x)\dx\ (2.7) 

Jo J-oo Jo 

poo poo p — ^ — 

/ / v t MwH(t -t-a 1 x)H(t-a 2 x)dtdx 4 / ^ \M(x)\ dx, (2.8) 

Jo J-oo Jo 

with the constant determined only by the upper bounds on |v|, |w|, | Vf| , w^| on the region {(x,t) : 
< <r 2 x < t < t — a\x\. 

For future use we note that since £5U = (5A)U X + (5B)U we observe that 

£(u - u) = (SA)u x + (5B)u, (2.9) 
£(u - u) = (5A)u x + (5B)u. (2.10) 

Also, from the construction of /(•) and £(•) we know that there are C 1 vector functions cf>(x,t) and 
<f)(x,t) so that 

(u + I * u)(x, t) = 6(t - nx)e 2 + cf)(x, t)H(t - px), (2.11) 
(u + 7*u)(z,i) = 5(t - nx)e 2 + 4>{x,t)H{t - fix) (2.12) 
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and the C 1 norms of and cf> on appropriate finite regions are bounded by A, fi, X and K. For 
future use we note that since u is given by we may conclude that 

u x (x,t) = -XS'(t - Xx)e 1 +fS(t - Ax) + (g- f)5(t - fix) +f x H(t - Ax) + (g x -i x )H(t - fix). 

(2.13) 

Below all constants are determined only by A, /i, X and K. 



Case I 

From (II. lip we have 

f3(x, t) = 5(r -t- Ax)ei + f (x, r - t)H(r - t - Ax) + (g - f)(x, r - t)H(r - t - fix). 

From (|2.9p we note that Col = (5A)u x + (5B)u. Now u, u x are given by (|2.ip and (|2.13p . so some 
important contributions to the LHS of (|2.3p from some singular terms in (3 T Ca are 

pOO pOO 

/ / ef (SA)ei S'(t - Ax) 6(t - t - Ax) dt dx = 0, 

JO J-oo 

(&B)ei<S(r - t - Xx)6(t - Ax) dt dx = (5B)n(x/(r))/(2A). 



f'OG /*00 

/ 



/0 J-oo 



All other terms on the LHS of (|2.3p may be estimated using (|2.5p - (|2.8p . Hence, using (5A)(0) = 0, 
we have 

\(5B) u (x f (T))\ 4 |RHSof &®\ + \(5A)(x f (r))\+ H * \{SA)(x)\ + \(SB)(x)\ dx 



o 



^ |RHS of ((23D|+ / ^ ] \(5A)'(x)\ + \(5B)(x)\dx. (2.14) 



/o 



Case II 

From ([2J)D, (pnUJ) we note that Ca = {8 A) (u + / * u) a . + (£B)(u + Z * u) so 



Ca(x, t) = (5A) y—fi5'(t — fix)e 2 — fi4>i x , t)5(t — fix) + (f> x (x, t)H(t — fix)j (2.15) 
+ (5B) (5(t - fix)e 2 + 4>(x, t)H{t - fix)) (2.16) 



and 

j3(x, t) = 5(t — t — fix)e 2 + 4>{ x i T ~ t)H( T — t — fix). 
Some important contributions to the LHS of (12. 3D from some singular terms in f3 T Co. are 

poo poo 

-fi I / ef (5A)e 2 8'(t - fix)6(r - t - fix) = 0, 

JO J-oo 

poo poo -| 

/ / 4(5B)e 2 5{t-fix)5(T-t-fix) = —{5B) 22 (x s (T)). 

Jo J-oo 2 V 
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All other terms on the LHS of (|2.3p may be estimated using (|2.5p - ([2.8p . Hence, as before, we have 

\(5B) 22 (x s (t))\ 4 |RHSof d23])| + r iT) \(5B)(x)\ + \(5A)'(x)\dx. (2.17) 

Jo 

Case III 

From ()2.11|) we have 

(3(x, t) = 5(t — t — fix)e 2 + <p(x, t — t)H(r — t — /ix) 

and from (|2,9p we have Ca = (5A)u x + (5B)u. Using 5,5^ given by (|2.ip . (|2.13p . some important 
contributions to the LHS of (|2.3p from some singular terms in fl T £a are 

roo /'OO ^ 

-X / e T(5A)(x)e 1 5(T-t- f ix)5'(t-\x) = - (<L4) 21 (x m (r)), 
/ / e 2 r (5B)e 1 5(t-\x)5(T-t-(ix) = -—(5B) 2 i(x m (T)). 

Jo J-oo A + H 

All other terms on the LHS of (12. 3p may be estimated using (|2.5p -( |2~?8j) . Hence, as before, we have 

\ fX m (T) 

\((SB) 21 -— — (8A)' 21 )(x m (T))\ 4 |RHS of <&M\ + / |(&B)(x)| + |(<L4)'(z)| dz. (2.18) 
A + Jo 

Case IV 

From (|1.11|) we have 

(3(x, t) =5{t -t- Ax)ei + f (x, r - t)H(r - t - Xx) + (g - f)(x, r - i)iT(r - t - fix) 

and Lex is the same as in Case II and is given by (|2.16p . So some important contributions to the 
LHS of (|2,3p from some singular terms in f3 T £ot are 

[ eJ{6A)(x)e 2 5{T - t - Xx)5'(t - fix) = ~^ (<L4)' 12 (x m (r)), 
Jo J-oo [a + M 2 

f'OQ roc -i 

/ / eJ(5B)e 2 5(t-^x)5(r-t-Xx) = -—(5B) 12 (x m (T)). 
Jo J-oo X + L 1 

All other terms on the LHS of (|2.3p may be estimated using (|2.5p - (j2.8p . Hence, as before, we have 

\((6B) 12 - -^—(5A)' l2 (x m (r))\ 4 |RHS of + / \(5B)(x)\ + \(5A)'(x)\dx. (2.19) 
A + pt Jo 

Fix an x in [0, XX/ fj] and define t s (x) = 2fix, t m (x) = (A + n)x, tf(x) = 2Xx to be the two-way 
travel time to probe a distance x at slow, mixed or fast speeds respectively. Then (I2.14p . (|2.17p . 
(fXTgj) . (J2HS1), together with (|2"1| may be combined into 

\((X + fi)(5B) - (5A)'VC)(x)\ 4 \(5U) x (0,t s (x))\ + \(5U) x (0,t m (x))\ + \(5U) x (0,t f (x))\ 

+ [ tsiX \(5U) X (0, t)\ dt+ [ X \(5A)'(z)\ + \(5B)(z)\dz 



QED 



^ max (5U) X (0,-) + / \(6A)'(z)\ + \(6B)(z)\ dz. 
[Q,ts(x)\ Jo 
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3 Proof of Theorem [2] 



Below all convolutions will be convolutions in the time variable only. Because of the ideas discussed 
in the introduction it is enough to prove Theorem [2] for the special case when diag(A) = - we will 
assume that for the rest of the proof. 

We must find an l(t) supported in [0, oo) so that v(x,t) := u(x,t) + l(t) * u(x,t) is zero on 
< Ax < t < fix. From ([l.lip . (|1.12|) we see that the most singular term in v(x,i) is 5(t — fix)e2 
but this has no impact in the region < Ax < t < fix. So, for the rest of the proof we will 
identify v(x, t) with v(x, t) — S(t — fix)e 2 over the region < Ax < t < fix. Now over the region 
< Ax < t < fix one may observe that 

v(x, t) = f (x, t) + l(t) * (5(t - Ax)ei + f (x, t)H{t - Ax)) 

rt— Xx 

= f(x,t) + l(t-Xx)e 1 + / l(s)¥(x,t - s)ds. (3.1) 
Hence we have to find an l(t) so that 

s) ds + fi(x,t) =0, 0<Xx<t<fix, (3.2) 
s) ds + f 2 (x,t) =0, 0<Ax<t</ix. (3.3) 
Fix a t > 0; then rewriting (|3.2p for points on the line t + Ax = r, we seek a function L(-) so that 

,-t-2\x 

L(t-2Xx)+ L(s)fi(x,T-Xx-s)ds + fi(x,r-Xx) = 0, x m (r) < x < x/(r). (3.4) 

Jo 

The Volterra equation (|3.4p has a unique solution L(-) in C[0, r(fi — X)/{fi + A)]. Since f and f 
are in C 2 and L is continuous, (j3.4j) implies that L e C 1 [0,r(fi — X)/{fi + A)] which again by (|3.4p 
implies that L G C 2 [0, r(/i — X)/{fi + A)]. The L(-) constructed depends on r but we have to find 
an L(-) independent of r - except for the domain of L which will depend on r. Moreover this L(-) 
must also satisfy (|3.3p . Both these goals will be achieved if we can show that vi(x, r — Ax) = for 
Xm(r) < x < Xf(r) implies that v(x,i) = for < Ax < t < fix, t + Ax < r; see Figure [3j Note 
that the supremum of L(-) on [0, t(/i — X/{fi + A))] is bounded above by a function of the supremum 
of /i(x, t) and /i(x, t) on the region 0<Ax<t<r — Ax. Hence, by Theorem [3l the supremum of 
L(-) on [0, r(fi — X)/(fi + A)] is bounded by a function of the supremum of A(-),B(-) on [0,xj(r)]. 

Below we use L{t) instead of l(t). From (|3.4p and (jl.2ip we observe that 

L(0) = -/ 1 (r/A,r) = ^=^oi a (0) 

so using ([33]) and (TOTD . (fl~22|) . we see that 

v(x, Ax) = f (x, Ax) + L(0)ei = 0. 



t—\x 



l(t - Ax) + / l(s)fi(x,t 



t~\x 



l(s)f 2 (x,t 
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t 




Figure 3: Interval of dependence for I 
Further, from (|3.ip . over < Ax < t < fix we have 

/•t—Xso 

v t (x, t) = f t (x, t) + L'(t - Xx)e 1 + L{t - Ax)f (x, Ax) + / L(s)f t (x, t- s)ds 

Jo 

so using p.!6p . (|1.22p we have 

vzt(x, Ax) = fzt(x, Ax) + L(0)/ 2 (x, Ax) 

= 777~~ ^-^ fl i2 (0)a 2 i(x) - ^ 9 ai 2 (0)-7p-^ — ,021 (») = 0. 
(A z — /i^J z A z — /i A z — /i z 

Hence from (jl.lip . (|1.12j) we see that 

Cvu — v xx — Av x — Bv = 0, < Ax < t < fix, t + Ax < r, (3-5) 
t>i(x,Ax)=0, V2(x,Ax) = 0, t>2i(x,Ax)=0, < x < xj(r). (3.6) 

We have to show that if v±(x,t — Ax) = for x m (r) < x < x/(r), that is i>i = on i?5 then 
v(x, t) = on the region ORS. This will follow from some energy identities - the only complication 
being the two velocities. One could also do this by setting v\ = on the relevant part of t = r 
instead of the t + Ax = r but one would not obtain the optimal interval of dependence results. 



t 




Figure 4: Two speed energy estimates 
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Fix a slope k strictly between A and \i and choose an arbitrary a G [t — (X + k)xf,r — (X + k)x m ] 
(the line t = kx + a cuts t = t between 7? and S 1 ). In Figured! for certain a, A will lie on OS 
instead of OR - the calculations below are simpler in this case. We have the identities 

2{c 2 w tt - w xx )w t = (c 2 w 2 + w 2 x ) t - 2(w t w x ) x , (3.7) 
2(c 2 w u - w xx )w x = 2(c 2 w x w t ) t - (c 2 wf + w 2 x ) x . (3.8) 

Using (13.8|) with w = v±, c = A and (|3.6p and that v\ (and hence Xvu — v\ x ) is zero on RS, we obtain 

/ / 2(\ 2 v ltt - vi xx )v lx dx dt 
J JOABS 

= \ 2\ 2 v lx vu + k{\ 2 v 2 lt + v\ x ) dx- I 2\ 2 v lx v lt + X{X 2 v\ t + vf x ) dx 
Jab Jos 

+ / 2X 2 v lx v u + ^{X 2 v\ t + v\ x ) dx - I X(Xv u - vix) 2 dx 
JOA JBS 

X(Xv lt + v lx ) 2 + {k- \){X 2 v\ t + v\ x ) dx+ I X(Xv lt + v lx ) 2 + X)(X 2 v 2 t + v\ x ) dx. 

AB JOA 



Hence 

(k - A) / (X 2 v 2 t + v 2 x ) dx< [[ 2(X 2 v ltt - v lxx )v lx dx dt. (3.9) 
Jab JJoabs 

Next, use fj3.7j) with w = V2, c = fi and (|3.6|) : also construct positive a, b with ab = k and a < /i 
and 6 < 1. Then we have 



// 

J Jc 



2{[i v 2 tt - V2 XX )v2t dx dt 

OABS 

y?v\ t + v\ x + 2kv 2t V2x dx- fi 2 V2 t + v\ x + 2Xv 2x V2t dx 
AB JOS 

,2„,2 , „,2 , ...... .. , / ,,2, ,2 



+ / \x v 2t + v 2x + 2{iV2 X V2t dx + / \i v 2t + f 2x - 2Af 2 tV2x dx 
'OA JbS 

{av 2t + bv 2x ) 2 + (fJ 2 - a 2 )vl t + (1 - 6 2 )t4, 

AB 

+ / {nv 2t + v 2x ) 2 dx + \ (fi 2 - X 2 )v 2 t + (Xv 2t - v 2x ) 2 dx. 
JOA JBS 



(/u 2 - a 2 )vj t + (1 - 6 2 )t>L dx < 2(fi 2 v 2 tt - v 2xx )v 2 t dx dt. (3.10) 

AB J J OABS 



Hence 



Finally, for i = 1, 2 and any P on j4B, using (|3.6p we have 



^ 2 (P) = ^ 2 (Q) + 2 / 
Jc 



ViVu dt < v\ + % (it. 
QP ^QP 

Hence 

/ v 2 dx< v 2 + v 2 t dxdt. (3.11) 

Jap J J oabs 
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If we define 



E(a) := I 
Jab 

then (|3.9p . (|3.10p . (|3.1ip may be combined to show that 



vt| 2 + |v x | 2 + Ivl 2 dx 



E(a) 4 I |v| + |vaj| + \(aiv 2x + + bi 2 v 2 )vi x \ + \(a 21 v lx + b 21 V! + b 22 v 2 )v 2t \ dx dt 

I OABS 

< II l v ! 2 + l v ^! 2 + l v <! 2 dx dt 

OABS 



<4 / E(s) ds. 

JT-{\+k)Xj(T) 

Hence E(a) = on [r — (A + k)xf(r),T — (A + k)x m (r)] by Gronwall's inequality. 

QED 



4 Proof of Theorem [3] 

We first prove the uniqueness. If there are two solutions of (jl.ip . (|1.2j) . (jl.4p in C 2 ([0, oo), 2?'(R)) 
then their difference is also a solution in C 2 ([0, oo), P'(R)) of fTT]) . but with u(0,t) = 0, 

u(0,t). Convolving this difference with any compactly supported smooth function of t, we have a 
smooth solution of this homogeneous initial boundary value problem and hence it will be zero by 
standard energy estimates. Since the convolution was with an arbitrary function of t, the difference 
of the two solutions must be zero proving the uniqueness part of Theorem [3l 

If we can construct u and u in the forms (jl.lip and (|1.12p with f , g being of C 2 regularity then 
u and u will be in C 2 ([0, oo), X>'(R)). Now we construct expansions for u and u with the properties 
mentioned in Theorem 



4.1 Progressing wave expansion 

If c is a constant, h(x,t) an arbitrary function and s(-) a distribution then one may show that 

C{hs(t - cx)} = s{t - cx)Ch + s'{t - cx)T{£, c)h + d t {s'{t - car) {C - c 2 /} h) . (4.1) 
where the first order transport operator T is defined as 

T(£, c)h := { (C + c 2 I)d t + 2cld x + cA] h . (4.2) 

We seek u and u in the form given by (jl.lip . (|1.12p for some arbitrary C 2 functions f, g, f,g 
which we assume are defined for all x, t. Of course the value of f , g, f , g only on the relevant parts 
will be needed to determine u, u. 
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From the boundary condition (jl.4p we see that 

g(0,t) = 0, g(0,t) = 0. (4.3) 

We now determine the conditions determining f and g. Using (jl.lip and (|4.ip we have 

= Cu = -S(t- Xx)Be 1 + XS'(t - Xx)Ae 1 

+ H{t - \x)Cf + 8{t - Xx)T(C, A)f + d t (5(t - Xx)(C - A 2 /)f) 

+ H(t - fix)C(g - f) + 5(t - fix)T(C, M )(g - f) + d t (5(t - nx){C - M 2 /)(g - f)). 

Sell we insist that f and g satisfy (|1.13|) : further, on t = Xx we insist that 

(C - A 2 /)f = —XAe\ (4.4) 

T{C,X)i = Be l (4.5) 



and on t = fix we insist that 



Now (g3]), (HH) give 



(C-^ 2 /)(g-f) = (4.6) 
T(£,/x)(g-f) = 0. (4.7) 



f 2 (x,Xx) = - 9 a 2 i(x), fi(x,nx) = gi(x,fj,x); (4i 
X z — fi z 



further (|4.5p implies 

(2A 2 / lt + 2Xf lx + A012/2) (x, Xx) = b n {x) (4.9) 

((A 2 + f i 2 )f 2t + 2Xf 2x + Xa 21 h)(x, Xx) = b 21 (x); (4.10) 

and (14.7|) implies 

((A 2 + ^ 2 )( 5l - /x)i + 2 M ( 5l - A), + M a 12 (<72 - f 2 )){x,fjLx) = (4.11) 

(2a* 2 (<72 - / 2 )t + M92 ~ f 2 )x + mMsi - /i))(x,/xx) = 0. (4.12) 

From (|4.8p we have (51 — fi)(x,fj,x) = 0, hence the transport equation (|4.12p implies that (g 2 — 
f 2 )(x,LLx) is constant. But g(0,i) = and /2(0, 0) = Aa 2 i(0)/(A 2 - fi 2 ). Hence 

/ , w x Aa 2 i(0) 

(92 ~ f2){x,fiX) = ^ __ ^ ■ 

Since (g\ — fi)(x,fj,x) = 0, taking its derivative and using it in (|4.1ip . we may conclude that 

(91 - fl)t(x,nx) = ^ /i A ^ 2 a 2 i(0)Qi 2 (x). 



2 One may show that the conditions imposed below are not just sufficient but also necessary to have £u = and 
£u = 
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Finally, using (|4.8p in (|4.9p we obtain 



d A 2 
2A— Ax)) = &n(x) + ^ 2 _ ^ 2 Ql2 ( :c ) Q 2i( g )- 



Integrating this and using /i(0,0) = g\ (0,0) = 0, we obtain 



h{xAx) = kL bn{z)dz+ 2(/-\i)i 



ai 2 (z) a 2 \{z) dz. 



We now determine the conditions characterizing f and g. Using (I1.12|) and ()4.ip we have 

= Cu = — S(t — fix)Be 2 + fj,6'(t — fix)Ae 2 

+ H{t - Xx)Cl + S(t - Ax)T(£, A)f + d t (5(t - Xx)(C - A 2 /)f) 

+ H(t - nx)C{g - f ) + <5(t - /ux)T(£, M )(g - f ) + d t (5(t - fxx)(C - fj, 2 I)(g - f)). 

So we insist that f and g satisfy (j!.13|) : further on i = Ax we insist that 

(C - A 2 /)f = (4.13) 
T(£,A)f = (4.14) 



and on t = fix we insist that 



(C- M 2 /)(g-f) = -/Me 2 (4.15) 
T(£, M )(g-f) = Be 2 (4.16) 



From (I4.13|) and (|4.15j) . we obtain 



f 2 (x,\x) = 0, (g x - fi)(x,fix) = -2—^a 12 (x), (4.17) 

from (|4.14p we obtain 

(2A 2 / W + 2Xf lx + Aa 12 / 2 )(x, Ax) = (4.18) 

((A 2 + ^ 2 )/ 2t + 2Xf 2x + Aa 21 A)(x, Ax) = 0, (4.19) 

and from (|4.16p we obtain 

((A 2 + ^ 2 )(si - h)t + M9i - fi)x + MOi2(52 - h))(x,nx) = b 12 (x) (4.20) 

{2 fJ 2 (g 2 - f 2 ) t + 2fi(g 2 - f 2 ) x + M a 2 i(^ - /i))(x,^x) = b 22 (x). (4.21) 

Using (|4.17p in (|4.18j) we conclude that /i(x,Ax) is constant. Now g(0, 0) = 0, so from (|4.17p . 
A (0,0) = a 12 (0)^/(A 2 - ii 2 ). Hence 

/i(x,Ax) = ^ 2 _^ 2 ai 2 (0). 
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Next using (|4. 17[) in (|4.19j) . we conclude that 



f 2t (x,Xx) = - ^ a 2 i(x)f 1 (x, Ax) = jj2~^j2 a i2(0)a 2 i(x). 



Also, from (|4.21j) and (|4.17p we conclude that 



d [i 2 
2/j,— ((g 2 - f2)(x,fix)) = 622(2;) + ^2 _ ^ 2 a M x ) a 2l( x ) 

which (with initial conditions) implies that 

_ - M f x 1 f x 

(52 - h){x,nx) = —3 2T / a 12 (z) a 21 (z) dz + — b 22 (z)dz. 

1{A — fJ, ) J z,fl Jo 

Hence to prove Theorem O we have to show that the initial boundary value problems (|1.13|) - 
(|1.18p and (|1.19j) - ()1.25j) have C 2 solutions over appropriate regions. Since A € C 2 [0,oo) and 
B G C 1 [0,oo) the right hand sides of all the zeroth order boundary and characteristic conditions 
are C 2 functions and the right hand sides of the first order characteristic conditions are at least 
C 1 functions. Further the compatibility conditions area also satisfied at (0,0). Hence Theorem [3] 
follows from Proposition [J] in subsection 14.21 



4.2 The Characteristic Boundary Value Problem 

Pick a constant T > and define the upper and lower regions 

U T := {(x,t) : < fix < t < T}, 

L T := {(x,t) : < Ax < t < fix, t<T}, 

Dt = Ut U L t . 

Proposition 4. Suppose A,B are in C^T/A], p G C 2 [0, T], q G C 2 [0,T/^], r G C 2 [0,T/A] ; 
s\ G C 1 [0,T//i], s 2 G C 1 [0, T/X] and satisfy the compatibility condition at (0,0), that is p(0)— r(0) = 
q(0). Also suppose that F,G are C 1 on Lt,Ut respectively. Then the Goursat problem 

£f = F in L T , £g = G in U T , (4.22) 

with the boundary conditions 

g(0,t)=p(t) forte [0,T], (4.23) 

(g - f)(x,/xx) = q(x), (gi - fi) t (x,fix) = si(x), for < x < T/fi, (4.24) 

f(x, Ax) = r(x), f 2t (x, Ax) = s 2 (x), for < x < T/X, (4.25) 

has a unique solution with f, g in C 2 . Further 

ll f lie 2 + llsllc 2 ^ IIpIIc 2 + llqllc 2 + ll r llc 2 + llsllc 1 + IIFIIc 1 + IIGIIc 1 ( 4 - 26 ) 
with the constant determined only by \\A\\ c i, ||-£>||ci and X,/j,,T. 
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Proof of Proposition [4] 

The uniqueness follows from the analysis in the ^4 = 0,-B = 0,F = 0,G = case below. We 
only give an outline of the proof of the existence part, highlighting the parts of the proof which 
are not standard. First, we explicitly write the solution of (|4.22p - (|4,25p for the special case when 
A = 0,B = 0,F = 0,G = 0. Then we use this special solution to reduce the original problem to 
the case where p = 0,q = 0, r = 0,s = which we deal with using a Volterra equation approach. 
For such problems, existence in H 1 may be derived by extending functions by and appealing 
to standard results for IBVP on the region x > 0, t > 0. However, these results will not give us 
the higher regularity of f , g because the total solution does not have this higher regularity across 
t = fix. Hence one has to appeal to techniques specialized to the problem under consideration. 

(A = 0, B = 0, F = 0, G = case) 

In this situation, the equations decouple so the problem reduces to studying characteristic boundary 
value problems for the wave equation. The derivation of the formulas for /2,<72 is easy enough; for 
(x, t) with Xx < t < fix, f2 is determined by the values of ji and jn on t = Xx - see the triangle 
PMN in Figure [5j Hence, the transmission condition (|4.24|) gives us 52 on t = fix; then for any 
point P(x,t) with < fix < t we have 32 = 9i{Q) + 92(S) — 52 The expressions for f2(x,t) 
and g2(x,t) consist of the values of T2,q2,P2 at linear combinations of x,t and the integral of S2 
over an interval with end points which are linear combinations of x, t. Hence the C 2 regularity of 
fi-, 92 follows quickly from the regularity of p, q, r, s. 



The derivation of the formula for f\,g\ is not so clear cut because now s = Xy is a characteristic 
and hence the value of /1 on s = Xy alone is not enough to determine f\ on < Ax < t < fix. An 
implicit method is needed and the boundary condition on x = and the transmission condition on 
t = fix now play a role. One starts with f\ and g± as sums of unknown functions of t — Xx and 
t + Xx and the required boundary and transmission conditions lead to the determination of the 
unknown functions. We will not write the long expression for fi(x,t) and g\{x,t) which consist 
of the values of r%,pi, q\ at linear combinations of x,t and the integral of si over an interval with 
end points which are linear combinations of x, t. Hence the C 2 regularity of /i,<?i follows from the 
regularity of p, q, r, s. 

(General case) 



Let 4>(x,t) and i/)(x,t) be the C 2 solutions of (l4T22D - (|4T25|) for the A = 0, B = 0, F = 0, G = 




y 







Figure 5: Constructing /2, 52 when A = 0, B = 0, F = 0, G = 
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case. Then f — cf>, g — if) is the solution of (|4.22p - (|4.25[) except with p = 0, q = 0,r = 0, s = 
and F, G replaced by F + Acf> x + Bcf), G + Aip x + Btp which are still C 1 functions on Lt and Ut 
respectively. Since <fi(x,t) and tp(x,t) are C 2 , we need to prove Proposition [2] only for the case 
when p = q = r = s = 0. 

For functions i(x,t) and g(x,i) defined over the regions < Ax < t < fix and < fix < t re- 
spectively we define, over the region < Ax < t, the piecewise function 



{f,g}(M) 



f (x, t) if < Ax < t < fix, 
g(x, t) if < fix < t. 



The value of {f , g} on t = fix is ambiguous and is to be understood to be the one sided limit. For 
any vector function H(x, t) on the region < Ax < t, we define a vector function T(H)(x, t) on the 
region < Ax < t as (see Figure [6]) 



J(H)(x,i) 



J JPQRS J JPLMN 



IPQRS 

where PQRS has sides parallel to s — Xy = or s + Ay = and PLMN has three sides parallel to 
s — fiy = or s + fiy = 0. Note that PLMN will change into a triangle PMN if Ax < t < fix. 







P(x,t) 










Q 










s+|J.y=0 


L 








s=Xy 


s+\y=0 








R 





Figure 6: Solution of inhomogeneous wave equation 



Now 

"t+Xx rt—Xx 



i pt+Ax rt—AX 

H 1 = — / H 1 ((q-p)/2X,(q + p)/2)dpdq, (4.27) 

ZA Jt-Xx JO 



'PQRS 

and 

^ rt+fix pt—fix 
I PLMN 2/i Jt-fix y_( (U _A) g /( A (+A) 



2 rt+fix /-t—fix 

H 2 = — / H 2 ({q-p)/2fi,(q+p)/2)dpdq (4.28) 

z f 1 Jt-ax J -(u-X)q/(u+X) 



with the lower limit of the q integral being (fi + X)\t — fix\/(fi — A) in the triangular PMN case, 
that is when t < fix. 

If H is continuous on the regions < Ax < t < fix and < fix < t (but may have jumps across 
t = fix) then X(H)(x, t) is at least C 1 on each of those regions. Clearly T(H)(x, t) is continuous on 
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< Xx < t; further, the first component of Z(H)(x,i) is C 1 on < Ax < t because its derivatives 
in directions parallel to s = Xy and s + Xy = are the integrals of H\ on PS and PQ respectively 
and these vary continuously with P even across t = fix. Also, from Figure El X(H) is zero on x = 
and t = Xx. Also, the first order derivatives of the second component of 1(H) on t = Xx are zero 
because the derivatives of the second component, in < Ax < t < fix, in directions parallel t = fix 
and t + fix = c are integrals along PN and PM (note L is not present in this case) . Finally, if H is 
C 1 on the regions < Ax < t < fix and < fix < t then T(H)(x, t) is C 2 on each of those regions. 
Hence the parts of 1(H) on < Ax < t < fix and < fix < t are the unique solution of (|4.22p - 
(|4.25p when A = 0, B = and p = q = r = s = with {F, G} replaced by H. 

Hence the f , g we seek are the solutions of the Volterra like integral equation 

{f , g}(x, t) = T({F, G})(x, t) + I({Af x + Bi, Ag x + Bg})(x, t). (4.29) 

Fix a T > 0; we wish to solve (|4,29p on Df. For any T E (0, T], let B be the Banach space of 
piecewise vector functions {f , g} with f a C 1 function on Lt and g a C 1 function on Ut with the 
C 1 norms. Now for any H G B we have shown that 1(H) is also in B; further using (|4.27p . (|4.28p . 
one may show that for (x, t) E Dt 

|J(H)(x,t)|, |^T(H)(x,t)|, |^J(H)(x,t)| ^ (T + T 2 )||H|| S . 

Here the derivatives of T(H)(x,t) and H(x,t) are to be understood to be one sided at points on 
t = fix; the constant is determined only by A, fi. Hence we can define the map K, from B to B with 

/C({f , g}) = T({F, G})(x, t) + l({Af x + Bf , Ag x + Bg})(x, t). 

For arbitrary {f , g} and {f , g} in B, we may observe that 

||/C({f,g})-/C({f,g})|| B <C(T + T 2 )(|A| c( ^ ) + |i?| c( ^ ) )||{f-f,g-g}|| B 

with C determined by A and fi, implying K, is a contraction for T small enough. Hence K, has a 
fixed point in B and we have proved the existence of the unique solution of (I4.29P for T > small 
enough. 

Now suppose we have solved (I4.29P for T = T\ for some < T\ < T; so we have f * , g* on 
which solve (|4.29p . For any T E [Ti,T] we redefine B as before except that we require that f = f * 
on Lt x and g = g* on Ut x - Define K, as before; then because f, g agree with f*, g* respectively on 
-Dti an d satisfy (|4.29p for T = T\, one may see that 

\\]C({1 g}) - £({f, g})|| B < C(|T -T 1 \ + \T- T x \*){\A\ c{Dt) + \B\ c{Df) )\\{l - f, g - g}|| B 

because after the subtraction and the cancelation of the contribution to the integrals over the region 
Dt ± , the remaining integrals are over subregions of T\ < t < T and no line parallel to t = ±Ax 
or t = ± fix in this region will have a length exceeding a constant times T — T\\ the constant 
determined by A, fi. Hence as before, IC is a contraction for T — T\ small enough. The important 
point is that there is a positive lower bound on T — T\ for which /C is a contraction and this lower 
bound is dependent only on |^4|c(d^), \B\c{p^,)i \ P*> an d T and is independent of T\. So repeating 
this argument we can construct the solution of (|4.29p over Df. 
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The solution constructed, {f, g}, is C 1 on Lt and Ut- However for such {f, g} the right hand 
side of (1P9D is C 2 on L T , U T . Hence {f , g} is C 2 on L T and U T . 

QED 
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